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Abstract 

The dynamic sructure factor as the basic quantity describing the 
collective excitations in a fluid is considered. We consider the cases 
of neutral and Coulombic fluids. The classical method of moments is 
applied to construct the dynamic structure factor satisfying all known 
sum rules. An interpolational formula is found which expresses the 
dynamic characteristics of a classical or quantum fluid system in terms 
of its static correlation parameters. The analytical results based on 
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the theory of moments are compared with Molecular dynamics data 
for various model systems. 

1 Introduction 

In the past there has been considerable interest in the time dependence of 
correlation functions or equivalently of the frequency dependence of structure 
factors. These functions has been studied in neutral and Coulomb fluids both 
theoretically and by molecular-dynamic simulations 0. Under a neutral 
fluid we understand here a fluid of particles interacting via a short-ranged 
potential. That means the termin neutral fluids includes such nonneutral 
systems as dusty plasmas and charged colloidal suspensions || where the 
interaction between the charged particles of one subsystem is screened by the 
motion of particles from another subsystem. 

Several approaches are devoted to the study of dynamic properties of 
strongly interacting fluid systems. In a rather incomplete list we mention 
the approaches in Refs. ||, |5| based on the memory function formalism, the 
approaches based on the theory of moments || |7[], and the approach based on 
the quasilocalized charge approximation [|J. It is interesting to note that all 
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these approaches succeeded by exploiting the method of collective variables 
|| in various modifications. 

This paper gives a short overview of the application of the method of mo- 
ments to the determination of dynamic properties of coupled fluid systems. 
As the main quantity describing the dynamics of a systems we consider the 
dynamic structure factor. The dynamic structure factor may be measured in 
scattering experiments. The peaks in the dynamic structure factor determine 
the collective excitations of the system. There may propagate different col- 
lective excitations depending on the type of the system (neutral or Coulomb). 
Generally speaking in neutral fluids we deal with sound modes whereas the 
plasma mode is a finite frequency mode. The different behavior of the modes 
is connected with the different behavior of the interaction potential Fourier 
transform at small wavenumbers k. In a neutral fluid the Fourier transform 
is finite, in the Coulomb case it diverges for k — > 0. 



2 Dynamic properties of neutral fluids 

We consider a system of N particles of one species with masses m and in- 
teracting via a pair potential V(r). The Fourier transform of the interaction 
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potential satisfies the inequality V(k = 0) < oo. The Hamiltonian of the 
neutral fluid reads: 

" = E^5lV(*^), (!) 

1=1 l^J 

where Pi is the ith particle momentum. In what follows we will use a classical 
notation, though all the calculations are easily generalized to the quantum 
case. 

Define the particle density and its Fourier transform 

n(r, = "*(*)), n k (t) = Y,e lk - x >V (2) 

i i 

,the density-density correlation function 

g(r,t) = (n(r,t)n(0,0)) . (3) 

and the dynamic structure factor 

S(k, u) = — I" S^-^gir, t)dt dr (4) 
2-nn J-oo 

In order to construct the dynamic structure factor as a central function for 
the determination of the dynamic properties of the system it is useful to 
consider the frequency moments of the dynamic structure factor: 

M n (k) = j_^ n S(u,k) = - ^—n k (t)n_ k (0)J (5) 
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Due to the parity of the structure factor all moments with odd numbers are 
equal to zero. The zeroth and second moments read 

M (k) = S{k), (6) 

M 2 (k) = -k B T. (7) 
m 



where S(k) = (1/N) \ n ] e n_] e / is the static structure factor of the fluid. The 
fourth moment includes particle correlations and reads, 

M 4 (k) = 3k i (k B T) 2 /m 2 + Mi \k) 
MT\k) = ^V(k)k B T + ^- 2 (8) 
£ [S(k + q)- S(q)] (k ■ qfk B T V{q) . (9) 

The Nevanlinna formula of the classical theory of moments expresses the 
dynamic structure factor 

ciu ^-Ittt, E n+i(k,z) + q n (k,z)E n (k,z) 

in terms of a function q n = q n (k, z) analytic in the upper half-plane Imz > 
and having a positive imaginary part there lmq n (k,uj + irj) > 0, rj > 0, it 
also should satisfy the limiting condition: (q n (k, z)/z) — > as z — > oo within 
the sector 9 < arg(z) < n — 9. The polynomials D n (and E n ) can be found 
in terms of the first In moments as a result of the Schmidt orthogonalization 



procedure. The first orthogonal polynomials read || 

D x = z, D 2 = z 2 -u\, D 3 = z{z 2 - u 2 2 ) , (11) 
E x = M , E 2 = M z, E 3 = M (z 2 + uf - u 2 2 ) , (12) 

where ojf(k) = M 2 (k)/M (k) and u;f(fc) = M 4 (k)/M 2 (k). Consider first the 
approximation n = 1 leading to the correct frequency moments Mq and M 2 . 
Using the Nevanlinna formula Eq. (|10D we obtain (q x = g 1>r + iqis), 



s(k,u) = ^l (13) 

vr [u; 2 - wf(fc) + q 1>r (k,u)u] +qf i (k,u)u 2 



We have no phenomenological basis for the choice of that function q x (z) which 
would provide the exact expression for S(k, z). We mention that the physical 
meaning of the function h\{z) = —iqi(z) is that of a memory function since 



from Eq. (13) it follows that the inverse Fourier transform of the function 
C(k, z) = (l/iif) /f^j S(k, oS)/(z — uj) obeys the equation 



at 2 



+ u?C(M) + / dshAk^-s) { ' ,S) =0. (14) 
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A simple approximation is to put the function qi(z) equal to its static value 
qi(z) = qi(0) = ivik, ) and Eq. (|TJP simplifies to the equation of a damped 
oscillator with frequency u)\ and damping constant v. 

d 2 C(k,t) _ „ , „ .dC(k,t) , x 

g t2 ^ + ^iC(k, t) + u(k, )^ Ll = • (15) 




generalization to the quantum case (T = 0) reads u>o{k) = 2 mS(k) EH " 



Consider now the long-wavelength behavior k — ► 0. In this case the static 
structure factor S(k — > 0) = nksTux is determined by the compressibility 
k t = — (1/V) (dV /dP) T . Then the dispersion relation reads 



which differs from the familiar dispersion equation for the sound wave by 
the factor 7. For a model of independent oscillators: c p = c v and 7 = 1. 
Therefore the above approximation for the static structure factor based on 
the Nevanlinna equation with n = 1 represents the model of independent 
damped quasiparticles. 

To go beyond this approximation one has to choose the 3-moment ap- 
proximation n = 2 in the Nevanlinna hierarchy reproducing the moments 
Mq, M 2 and M 4 . Within this approximation and choosing q 2 (k,uj) = h(k) 
we obtain the following expression for the dynamic structure factor: 




(16) 



S(k,cu) 



S(k) h(k)ul(k) (ujl(k) - ujf(k)) 



(17) 



2 ■ 
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where h(k) has to be taken from the relation 

h(k) = (col - u*)/u(k) = (S(k)/K)((ul/ui - l)/S(k, 0)) (18) 

in order to satisfy the exact low freqency behavior S(k, 0). The value S(k, 0) 
may be taken from elastic scattering experiments, from another theory or it 
may be used as a fit parameter. 

Consider again the long wave-length limit k — > 0. Then the frequencies 

« 2 * 2 , ^=(^) (19) 

uj 2 2 (k)=v 2 k 2 , v 2 = -V(0) + 3^. (20) 

mm 

At small temperatures /c#T <C nV(0) we have u 2 = v 2 and we obtain the 
dynamic structure factor for a "classical" fluid at low temperature 

S(k,u) = ^ Y {5{u-ku) + 5{uj + ku)} , (21) 

representing undamped sound waves. The corresponding generalization to a 
quantum fluid reads 

S ^ = M.-eJwW) ^ - ^ + e "' k "' tBT ^ + *">} • (22> 
At zero temperature the system may only absorb energy and we obtain the 
simple equation for the dynamic structure factor. 

S(k,u) = —5(u-ku). (23) 
mu 



3 Dynamic properties of Coulomb fluids 

Consider a one component plasma (OCP) consisting of N particles with 
charges Ze and masses m interacting via the Coulomb potential V c (r) = 
Z 2 e 2 /r and embedded in a neutralizing homogeneous background. The clas- 
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sical OCP may be characterized by the coupling parameter T = a £ gT , a = 
(3V/4:irN) 1 ' 3 being the Wigner-Seitz radius. The quantum plasma has an ad- 
ditional parameter - the degeneration parameter 9 = a/ 2mA; #T /h 2 (37r 2 n) 2 / 3 . 
In what follows for the case of simplicity we concentrate on a classical plasma. 
For r < 1 we deal with an ideal (or Vlasov) plasma, for T ^> 1 the plasma is 
called a strongly coupled one. The Vlasov approximation takes into account 
only the mean field part of the interaction and the dispersion relation for the 
longitudinal plasmons is predicted as oj 2 {k) = u 2 ^1 + with the plasma 

frequency to 2 = ^ z lf n and the squared inverse Debye- length k 2 D = 47r ^y n ■ 
The Vlasov theory predicts a strong positive dispersion of the plasmons, i.e., 
cLu/dk > 0. However, in a coupled plasma the potential energy plays an im- 
portant role and the Vlasov approximation is not longer valid. To construct 
the dynamic structure factor for a coupled plasma consider the frequency 
moments of the dynamic structure factor S(k,u). The frequency moments 
formally coincides with that of a neutral plasma (Eqs. The only dif- 
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ference is that the interaction potential of the neutral fluid has to be replaced 
by that of the Coulomb system. The application of the Nevanlinna formula 
leads then to a corresponding hierarchy of approximations for the dynamic 
structure factor. If one is interested in the structure factor of a quantum 
system again Eqs. ( |10|) hold, if one replaces S(k,u) on the left hand side of 
the Eqs. (|T0|) by the loss function R(k,u = [(1 — exp(— (3%uj)) / (3Tilo\ S{k,uj). 
However, in the quantum case additional contributions to the zeroth and 
fourth frequency moment occur || [7j . 

Consider the 3-moment approximation Eq. fllTI). If one is interested in 
the investigation of the high-frequency collective excitation spectrum only 
it is sufficient to neglect the function h{k) since the damping (described by 
the function h) is small in strongly coupled plasmas. If one puts h(k) = 
Eq. (|Tj]) provides the expression of the dynamic structure factor for a 
strongly coupled plasma obtained within the QLC approach ||, if the thermal 
contributions may be neglected with respect to the correlation contributions. 
Within the simple approximation h(k) = the dynamic structure factor has 
S peaks at the frequencies uo c which in the classical case are determined by 
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the equation 



2/,\ ^4 2 



/l + 3 W + Jt E + (24) 



For — > the dispersion relation simplifies and we get 

with £7 C being the correlation energy density. Using the simple estimation 
k ^Tn = ~ 0.9r valid in the strong coupling regime one obtains the disper- 
sion relation u> 2 (k) = u 2 [1 + k 2 a 2 (— 0.08 + T -1 )] and one predicts a negative 
dispersion for T > 13. 

To study the dynamic structure factor one has to go beyond the simple 
approximation h = 0. To satisfy the low frequency behavior one may choose 
the approximation Eq. ( |I8|) . To check the quality of the predictions from 
our approximation molecular dynamic simulations have been performed for 
comparison J/J. The semiclassical simulations were performed to model a 
quantum gas of 250 electrons moving in a cubic box with periodic bound- 
ary conditions. The thermal equilibrium was established by a Monte Carlo 
procedure. A detailed description of the semiclassical model used in the 
simulations may be found elsewhere 0. In Figs. [3] and || we have plotted 
the loss function R(q,u) (q = ka) for various values of wavenumbers q for 
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the cases of strong (T = 10 ) and very strong coupling (T = 100 ) [0]. In 
both cases we obtain a sharp plasmon peak at small q values, with increas- 
ing wavenumber the plasmon peak widens. Almost no dispersion has been 
observed at V — 10. This is in good agreement with the above estimation 
for the critical value T = 13 separating regimes with positive dispersion from 
that with negative dispersion. For the case of very strong coupling T = 100 
we have found a strong negative dispersion. In Figs. |3| and f| we present the 
results of the MD data and compare them with our analytical approxima- 
tion Eqs. ( |TT| ) and (|1^). To calculate the parameters Ui(k) and uj 2 (k) we 
have used the static structure factor obtained from the HNC equations. The 
value S(k, 0) determining the parameter h(k) might be taken from the MD 
simulations. However, the dynamic structure factor at the zero frequency 
can be obtained with the necessary accurazy only from long time simula- 
tions. Alternatively we have choosen the value S(k, 0) to fit the model to the 
MD data. It should be mentioned that the value S(k, 0) mainly determines 
the width of the plasmon peak, the peak position is quite insensitive to the 
choice of the value S(k,0). From the figures it can be seen that there is a 
reasonable agreement between the MD data and the present approxiamtion 
based on the sum rules. The peak position is reproduced with high accuracy, 
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the agreement in the width of the peaks is less satisfactory. One concludes 
that the static approximation q2(k,u>) = ih(k) undersetimates the damping 
of the quasiparticles. 



4 Conclusions 

In this paper we have shown that the application of the classical theory of 
moments gives a satisfactory description of many properties of neutral and 
Coulomb fluids. The Nevanlinna formula generates approximate expressions 
for the dynamic structure factor in terms of their static correlations. The 
quality of the Nevanlinna expression mainly depends on the quality of the 
model used to calculate the static properties of the fluid. The presented 
results may be improved by a specification of the interpolation function 

In conclusion, the present approach has been also used to calculate the 



dynamic structure factor of two-dimensional electron gas [11|, of binary ionic 



mixtures [12] and of two-component plasmas [13]. It had been extended to 



magnetized plasmas [14]] and can be generalized to calculate partial dynamic 



structure factors. Here, the matrix form of the Nevanlinna formula becomes 



13 



helpful. 
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FIGURE CAPTIONS 

(Figure 1) The simulation data for the loss function R(q,u) versus fre- 
quency u/ujp for different wavevectors q = ka at T = 10 and 9 = 1. 

(Figure 2) Same as in Fig. 1 at T = 100 and 9 = 50. 

Figure 3 Comparison of the loss function R(q, uj) within the present sum 
rules approach (Eqs. ([L7| and (|18|) with S(k,u) replaced by R(k,u)) 
versus frequency uj/uo p with the corresponding MD data at V = 100 
and 9 = 50 for wavevector q = 0.619. 

Figure 4 Same as Fig|| at T = 100 and 9 = 50 for wavevector q = 1.856, 
01 • 
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Figure 1: The simulation data for the loss function R(q, uj) versus frequency 
u/ujp for different wavevectors q = ka at T = 10 and 9 = 1. 
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Figure 2: Same as in Fig. 1 at T = 100 and 9 = 50 
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Figure 3: Comparison of the loss function R(q,u>) within the present sum 
rules approach (Eqs. (|TT| and (p~8|) with S(k,u) replaced by R(k,u)) versus 
frequency uj/uj p with the corresponding MD data at T = 100 and 6 = 50 for 
wavevector q — 0.619. . 
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Figure 4: same as Fig.|3]; at T = 100 and 9 = 50 for wavevector q = 1.856, 
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